(N 
O 
O 
(N 



Stability of sub-surface oxygen at Rh(lll) 



to appear in Phys. Rev. B 



M.Veronica Ganduglia-Pirovancfl, Karsten Reuter, and Matthias Schefner 
Fritz- Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany 

(February 1, 2008) 

Using density-functional theory (DFT) we investigate the incorporation of oxygen directly be- 
low the Rh(lll) surface. We show that oxygen incorporation will only commence after nearly 
completion of a dense O adlayer (0 to t ~ 1-0 monolayer) with O in the fee on-surface sites. The 
experimentally suggested octahedral sub-surface site occupancy, inducing a site-switch of the on- 
surface species from fee to hep sites, is indeed found to be a rather low energy structure. Our results 
indicate that at even higher coverages oxygen incorporation is followed by oxygen agglomeration 
in two-dimensional sub-surface islands directly below the first metal layer. Inside these islands, 
the metastable hep/octahedral (on-surface/sub-surface) site combination will undergo a barrierless 
displacement, introducing a stacking fault of the first metal layer with respect to the underlying 
substrate and leading to a stable fcc/tetrahedral site occupation. We suggest that these elementary 
steps, namely, oxygen incorporation, aggregation into sub-surface islands and destabilization of the 
metal surface may be more general and precede the formation of a surface oxide at close-packed late 
transition metal surfaces. 
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PACS numbers: 81.65.Mq, 68.43. Be, 82.65.My 
I. INTRODUCTION 

Exposing a metal surface to oxygen can result in simple 
adsorbate covered surfaces, O sub-surface penetration or 
oxide formation, depending markedly on the oxygen par- 
tial pressure, substrate temperature, the crystallographic 
orientation of the surface, and the time of exposure. As 
the chemical activity of the metal surface in these dif- 
ferent states can vary significantly, knowledge of the in- 
teraction of oxygen with metal surfaces is critical for un- 
derstanding technologically important processes like e.g. 
oxidation catalysis. The formation of surface oxides on 
metal surfaces can in this respect both be beneficial, as 
well as detrimental. In particular, for applications such 
as the CO oxidation on transition metal based catalysts, 
oxide formation at the catalyst's surface in the reactive 
environment was primarily viewed as leading to an inac- 
tive surface oxide poisoning-tbje catalytic reaction, as e.g. 
in the case of Rh surfacesjjtil while recent experiments 
showed that the high catalytic activity of the Ru(0001) 
surface with respect to the CO oxidation relates in fact to 
the existence of Ru02(U(lj oxide patches that form un- 
der oxidizing conditions UU Thus, apparently, unreactive 
surface oxides are more likely to form on Rh than on Ru 
surfaces. These findings call for an atomistic understand- 
ing of the oxidation of these transition-metal surfaces and 
of the structure of their oxide surfaces. 

There is general agreement that the mechanistic steps 
leading to oxide formation involve dissociative oxygen 
chemisorption on the metal surface, lattice penetration 
of atomic oxygen, and the crystallization and growth of 
the stoichiometric oxide phases. The sequence of these 
events may be complex, and rather than successively, 
they may occur simultaneously, depending on the tem- 
perature and pressure. The first of these steps, i.e. the 
chemisorption of oxygen on Rh single crystal surfaces, 



has been the subject, of intensive experimental and the- 
oretical research,D~li3 but only few experimental stud- 
ies have-addressed the ensuing formation of sub-surface 
oxygen JL1 i3. In one of these studies it could recently 
be showr£3 that exposure of the Rh(lll) surface to O2 
at moderatly elevated temperatures (~ 470 K) leads to 
the formation of sub-surface oxygen species, and there 
is evidence that at least ~ 0.9 monolayer (ML) of O is 
adsorbed on the surface before sub-surface sites are occu- 
pied. Interestingly, the analysis of the X-ray photoclcc- 
tron diffraction (XPD) data of the corresponding study 
suggested that the sub-surface species (~ 0.1 ML) oc- 
cupies octahedral sites between the first and second Rh 
outermost layers, which lie just underneath the fee on- 
surface adsorption sites, while the neighboring on-surface 
oxygen has switched from its normal fee to the hep sites. 

On theory side, ab initio studies of O sub-surface 
species-, on. transition metal surfaces are even more 
scarceliaErj and as yet lacking for Rh surfaces. The 
present theoretical study therefore specifically examines 
the incorporation of O into the close-packed Rh(lll) 
surface using density-functional theory. Extending pre- 
ceding rW°fk focusing on ordered oxygen phases on this 
surfacdaM we now discuss the stability of oxygen be- 
low the surface, and address questions concerning the 
minimum oxygen coverage at the onset of oxygen pene- 
tration (section III. A), the stability of various available 
sub-surface sites (section III.B), as well as the effect of a 
continued oxygen incorporation with increasing coverage 
up to 2.0 ML (section III.C). The present results for sub- 
surface oxygen on Rh(lll) are discussed and compeared, 
with corresponding results for the Ru(0001) surfacal§E£l 
as a first step towards the understanding of the elemen- 
tary steps governing the onset of surface oxide formation. 
Questions remain as to the nature of the distinct chemi- 
cal activity of oxidized Ru(0001) and Rh(lll) surfaces. 
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II. CALCULATIONAL DETAILS 

All calculations have been performed using density- 
functional theory (DFT) and the generalised gradi- 
ent approximation (GGA) of Perdew et alE3 for the 
exchange-correlation functional as implemented in the 
full-pot&H|ti|aJ linear augmented plane wave method (FP- 
LAPW)ErM The Rh(lll) surface is modeled in the su- 
percell approach, employing a 7-layer (111) Rh slab with 
a vacuum region corresponding to 6 interlayer spacings 
(w 13 A). Oxygen atoms are adsorbed on both sides of 
the slab. Their positions as well as those of all Rh atoms 
in the two outermost substrate layers are allowed to relax 
while the central three layers of the slab are fixed in their 
calculated bulk positions. In a preceding publication we 
have detailed the geometrical properties and stability of 
ordered adlayers of O on Rh(lll)B. Exactly the same cal- 
culational setup is used here, such that the results of our 
prior study can be used for comparison. Further tech- 
nical details of the calculations can correspondingly be 
found in Ref. g 

We address the stability of 0/Rh(lll) structures with 
respect to adsorption of O2 by calculating the average 
adsorption energy per O adatom, 

E ad (9 tot )=E b (6 tot )-D/2, (1) 

where D is the dissociation energy of the O2 molecule and 
Eb(Otot) the average binding energy per oxygen atom as 
a function of the total coverage (i.e., on- + sub-surface) 
#tot- In turn, i?b(0tot) is defined as 

E h (0tot) = I S O/Rh(in)(0tot) ~ ( 2 ) 

(^Rh(lll) + ^totE Q ) \ , 

where iVtot is the total number of oxygen atoms in the 
unit-cell, and £o/Rh(m)> S Rh(m)> and £ o om are the 
total energies of the 0/Rh(lll) adsorbate system, of the 
clean Rh(lll) surface, and of the free oxygen atom, re- 
spectively. A positive value of the average adsorption 
energy indicates that the dissociative adsorption of O2 is 
exothermic. That is, the binding energy per O adatom 
on Rh(lll) is larger than that which the O atoms have 
in 02(gas), i.e., D/2. The total energies of the adsorbate 
system and clean surface are calculated using the same 
supercell. Details on the calculations of the isolated O 
atom and the free O2 molecule are given in Ref. ||. 

To test the accuracy of the calculated binding ener- 
gies, Eb(9 to t), cm numerical approximations due to the 
finite FP-LAPW basis set and the finite slab and vac- 
uum thicknesses in the supercell approach, selected cal- 
culations were repeated with higher accuracy. The self- 
consistent calculations of E b were routinely conducted 
for a 7-layer slab, a 16 Ry plane-wave cutoff in the inter- 
stitial region, and a (12 x 12 x 1) Monkhorst-Pack grid 
for thft,(l x 1) unit-cell with 19 k-points in the irreducible 
wedgaa. Changing to denser k-meshes up to a (18 x 18 x 1) 



grid with 37 k-points in the irreducible wedge resulted in 
negligible variations of E b within ~ 10 meV/O atom. 
Similarly, extending the Rh(lll) slab from 7 to 9 and 11 
layers led only to changes in the calculated binding en- 
ergies up to ~ 10 meV/O atom. The only really notable 
effect on the computed E^ is caused by variations of the 
finite plane-wave cutoff in the interstitial region: Increas- 
ing this cutoff from 16 Ry to 24 Ry decreases the binding 
energies by w 100 meV/O atom. Yet, this decrease can 
be largely attributed to an improved description of both 
the free O atom and the chemisorbed species and thus 
affects all structures that contain the same amount of O 
alike. Consequently, the calculation of relative stabili- 
ties of different phases, i.e. the difference between two 
binding energies, has a significantly smaller error than 
this variation in the absolute value of E^. Combining all 
these tests, we give a conservative estimate of this latter 
error of ±30 meV/O atom, which in turn does not affect 
any of the conclusions made in this work. 

III. STABILITY OF SUB-SURFACE O 
A. Minimum coverage for O incorporation 

The initial oxidation of the Rh(lll) surface proceeds 
via the dissociative chemisorption of oxygen in the three- 
fold fee hollow sites of the basal surface.LTLl O2 exposure 
under UHV conditions leads to two well ordered super- 
structures, namely a (2 x 2)-0 and a (2 x l)-0 phase at 
coverages of 0.25 ML and 0.5 ML, respectively. The lat- 
ter coverage was for,a long time believed to correspond|-tp. 
a saturated surface.Q Yet, after theoretical predictionJj'El 
a (1 x l)-0 phase (0 to t = 1.0 ML) could be stabilized 
experimentally under UHV conditions using atomic O as 
oxidant, thus demonstrating that the apparent saturation, 
results from a kinetic energy barrier to O2 dissociation.EZI 
This barrier can also be overcome by eahajaeed O2 dosage 
at moderately elevated temperaturesEJOcj or by eptap-. 
sure to more oxidizing carrier gases like e.g. N02»'tLZl 
In all cases, sub-surface O incorporation or bulk diffu- 
sion was reported for higher exposures at temperatures 
above ~ 400 K, and as already mentioned a recent exper- 
imental study brought evidence that at least ~ 0.9 ML 
of oxygen is adsorbed pa the surface, before sub-surface 
sites become occupied.E£l 

To address this initial incorporation theoretically we 
study the stability and properties of sub-surface oxygen 
species as function of the total coverage, to t (0.25 < 
#tot < 1.0 ML). We employ (2 x 2) unit-cells and calcu- 
late the average adsorption energy of fully relaxed struc- 
tures containing from N to t = 1 (0tot = 0.25 ML) up to 
Y tot = 4 (0 tot = 1 ML) oxygen atoms (see Eq. (§)). The 
occupation of sub-surface sites will commence, when a 
structure with one of these O-atoms located below the 
surface becomes energetically more favorable than the 
most stable one with all oxygens on the surface at the 
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hcp/tetral fee /tet rail (a) fcc/tetral I (b) hcp/tetrall (a) hcp/tetrall (b) 

FIG. 1. Top view of all possible on-surface/sub-surface site 
combinations at f5 to t = 1.0 ML with one oxygen atom located 
below the surface (see text for the explanation of the differ- 
ent sites). Rh = big spheres (white = surface layer, grey 
= 2nd layer), O = small spheres (black = on-surface, grey = 
sub-surface). Oxygens in tetra-II sites below the first layer Rh 
atoms are invisible in this plot and are schematically shown 
as small white circles. Symmetry inequivalent occupation of 
the same kind of on- and sub-surface sites are denoted with 
(a) and (b), respectively. 



same total coverage. That is, we compare the stability of 
(2 x 2) - ((iV tot - l)O on + O sub )/Rh(lll) mixed phases 
relative to the (2x2) — (Ar tot O ii) /Rh(lll) structure with 
O on atoms in fee hollow sites for iV to t = 1, . . . 4. 

Between the first and second outermost metal lay- 
ers, there are three different high-symmetry interstitial 
sites available for O incorporation. The octahedral site 
(henceforth octa) lies just underneath the fee on-surface 
site, and one tetrahedral site (tetra-I) lies below the hep 
on-surface site. A second tetrahedral site (tetra-II) is lo- 
cated directly below a first layer metal atom. Considering 
both three-fold hollow sites (fee and hep) as possible ad- 
sorption sites on the surface, leads then to a many-fold of 
possible structural combinations of how to place the iV t ot 
oxygen atoms into the unit-cell, in particular as at most 
coverages several symmetry inequivalent possibilities for 
the same on- and sub-surface site combination exist. 

As will be shown below, the stability of structures in- 
volving on- and sub-surface sites becomes only compa- 
rable to that of the pure chemisorbed on-surface phase 
at #tot = 1.0 ML. Thus, only at that coverage we will 
address the complete set of possible structural combina- 
tions, shown in Fig. [j], while at the lower coverages we 
only exemplify the stability trends by computing three 
likely combinations that will become relevant at a later 
stage of our discussion. Namely, these are geometries 
where on-surface oxygen is in fee sites and the sub-surface 
oxygen in either the octa or the tetra-I sites, fcc/octa and 
fcc/tetra-I respectively. Thirdly, we included the exper- 
imentally suggested possibilityta that oxygen in octahe- 
dral sub-surface sites induces a site-switch of the nearby 
on-surface oxygens from fee to hep, viz hcp/octa. The 
thus resulting set of considered structures is shown in 
Fig. | for to t = 0.5 ML and tot = 0.75 ML. At the 
latter coverage we also tested the possibility of a simul- 





fcc/octa hcp/octa (a) hcp/octa (b) fcc/tetral (a) 
(2 x 2)-(20 on + O sub )/Rh(lll), 6 tot = 0.75 ML 

^v/. :* • v • a *. 

fcc/octa fee hcp/octa hcp/octa (a) hcp/octa (b) fcc/tetral (a) 

FIG. 2. Top view of selected on-surface/sub-surface site 
combinations at f?tot = 0.5 ML and f5 to t = 0.75 ML with one 
oxygen atom located below the surface (see text). Rh = big 
spheres (white = surface layer, grey = 2nd layer), O = small 
spheres (black = on-surface, grey = sub-surface). 

TABLE I. Average adsorption energies (in eV/O atom) 
of Rh(lll) geometries containing O in on-surface and/or 
sub-surface sites. The label on/sub indicates the site type 
(fee, hep for on-surface; octa/tetra-I/tetra-II for sub-surface, 
see text). Calculated values for the most stable on-surface ad- 
layer with O occupying fee sites (fee/ — ) are from Ref. For 
0tot < 1.25 ML, 0.25 ML of O is contained below the surface, 
for t ot = 2.0 ML, 1.0 ML of O is contained below the surface. 
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taneous occupation of hep and fee on-surface sites and 
octahedral sub-surface sites, leading to the fee hcp/octa 
structure. Figs. |l] and || also show the considered sym- 
metry inequivalent structures with the same kind of on- 
and sub-surface sites occupation, e.g. hcp/octa (a) and 
hcp/octa (b). 

The calculated average adsorption energies of all in- 
vestigated structures are compiled in Table |, while Fig. 
^ additionally visualizes the trends for three selected 
on-surface/sub-surface combinations as described above. 
Concentrating first on the lowest tested coverage, 9 to t — 
0.25 ML, we find that the occupation of just sub-surface 
sites without any on-surface oxygen is not even exother- 
mic and by « 3eV/0 atom less favorable than adsorp- 
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FIG. 3. Average adsorption energies (in eV/O atom) of se- 
lected Rh(lll) geometries containing O in on-surface and/or 
sub-surface sites, see text and Table fl. 



face can only occur at just about the completion of the 
full monolayer coverage on the surface, i.e. 6 to t ~ 1.0 ML. 

The calculated energy difference of « 0.3eV/O atom 
between the pure fee and the mixed phases for cover- 
ages 0tot ~ 1-0 ML indicates that sub-surface O pene- 
tration into Rh(lll) is an energetically activated process 
and that a small, but finite concentration of sub-surface 
oxygen can be expected already at on-surface coverages 
slightly below 1.0 ML at elevated temperatures. Wider 
et al. have found that extended exposure of the Rh(lll) 
surface to oxygen at 470 K led to the formation of a sub- 
surface species, though-the on-surface coverage remained 
slightly below 1.0ML.EJ3 In comparison, on Ru(0001) we 
find the energy difference between the pure on-surface 
and the pixed phases at #tot = 1.0 ML to be » 0.8 eV/ 
O atom,EZl that is about a factor of three larger than 
for Rh(lll). Correspondingly, oxygen penetration has 
at the prior surface only been reported to occur after th£ 
(1 x l)-0 phase on the surface has been completed. t3E9 



B. Site preferences and site-switch 



tion into the on-surface fee hollow sites. We have recently 
shown that this is generally the case for the closed-packed 
late Ad transition metal surfaces, which is largely largely 
due to the significant local expansion of the metal lat- 
tice induced by occupation of sub-surface sites.EEl The 
corresponding cost of distorting the metal lattice and 
breaking metal bonds renders sub-surface sites initially 
always less stable than on-surface chemisorption. Yet, 
upon increasing the on-surface coverage, repulsive inter- 
actions between the adsorbates drive the adsorption en- 
ergy down, cf. Fig. [jj, until eventually occupation of 
sub-surface sites might become more favorable compared 
to a continued filling of on-surface sites. 

Turning therefore to higher coverages up to 1.0 ML, 
where we continuously increase the number of adsorbed 
on-surface oxygen atoms per unit cell, we find the aver- 
age adsorption energies of all mixed structures now to 
be positive, i.e. they should be able to form. Still, in 
the whole sub- monolayer regime the values of 2? a d(0tot) 
for the mixed phases are significantly lower compared to 
that of the pure on-surface adsorption at the same total 
coverage, reflecting the above described fact, that the oc- 
cupation of sub-surface sites is energetically considerably 
less favorable. Only at 8 to t — 1-0 ML does the adsorp- 
tion energy of the selected mixed phases approach that 
of the pure on-surface fee phase to within w 0.3eV/O 
atom (cf. Fig. ||). To make sure that no other struc- 
tural on-surface/sub-surface combination would be even 
slightly more stable and then eventually more favorable 
than the pure on-surface phase, we tested all of the possi- 
ble structures shown in Fig. Iflat this particular coverage. 
As can be seen from Table [j] neither of these configura- 
tions leads to a more favorable binding than the pure fee 
phase. Therefore, O incorporation into the Rh(lll) sur- 



Trying to understand the oxygen site preferences in the 
various on-surface/sub-surface phases at coverages up to 
1.0 ML, we note first that the low stability of a number 
of structural combinations can be understood in terms of 
electrostatic repulsion between the oxygens, i.e. when- 
ever the electronegative oxygens come too close to each 
other. This holds for mixed fee hep combinations (cf. 
Table |) , as well as for the less stable of the two symme- 
try inequivalent possibilities of the same on-surface/sub- 
surface site combination (rendering e.g. the hep/octa (b) 
and fec/tetra (b) geometries unfavorable, cf. Figs. |l| and 
| and Table J). 

Moreover, we find the stability of all other combina- 
tions at #tot = 1-0 ML to be rather similar, i.e. within 
w 0.2eV/O atom, cf. Table I. This suggests that at ele- 
vated temperatures oxygen incorporation could start ini- 
tially in all available sites, and that kinetic factors (like 
penetration barriers) more than energetic factors deter- 
mine which of the possible metastable sites get populated 
first in an experiment with controlled oxygen dosage. As 
noted in the introduction, a recent analysis of XPD data 
indicated the presence of a sub-surface oxygen species 
(~ 0.1 ML) in octahedral sites between the first and sec- 
ond metal layer at an almost oxygen-covered (~ 0.9 ML) 
Rh(lll) surface after 10 5 Langmuir O2 exposure at 470 
K. The analysis of the data yields to the suggested simul- 
taneous occupation of neighboring on-surface hep sites, 
which is different from the otherwise preferred fee sites 
on Rh(lll) before oxygen incorporation.tj 

This intriguing finding is in line with our calculations, 
which indeed indicate, that if sub-surface oxygen was to 
occupy octahedral sites, then the site-switch combina- 
tion hep/octa would be considerably more favorable com- 
pared to the fec/octa combination with the on-surface 
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FIG. 4. Calculated work function change for the fec/octa 
and hep/octa (a) structures at f?tot = 1.0 ML (see Fig. Cal- 
culated values for the most stable adlayers with O occupying 
only on-surface fee sites are from Ref. pi 
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FIG. 5. Calculated initial-state on-surface O Is 

core-level shifts for the fec/octa and hep/octa (a) phases at 
dtot = 1-0 relative to the O Is level for the adatoms of the 
0-(l x 1)/Rh(lll) structure at 6 on = 1.0 ML, with O occupy- 
ing on-surface fee sites (fee/ — ). The label on/sub indicates 
the on- and sub-surface adsorption sites, respectively. There 
are two symmetry inequivalent O adatoms in the hep/octa 
phase, which is why two levels are shown for this structure. 

oxygens in their normal adsorption sites, cf. Fig. ^. 

In seeking for a qualitative explanation for the stabil- 
ity of the hep on-surface sites upon occupation of octa- 
hedral sites, and consistent with the discussion on Iha 
preferred fee site adsorption begun in earlier papersJj'EHI 
we examine the calculated changes in the work-function, 
A$, as well as the difference in O Is binding energy of 
on-surface O-atoms, AOi s , of the mixed fec/octa and 
hep/octa phases at # to t = 1-0 ML relative to the values 
for the stable (1 x l)-0/Rh(lll) structure, with O just 
occupying fee sites. The calculated A$ and AOi s values 
are shown in Figs. ^ and respectively. Initially, the 
work function rises as a function of the coverage, while 
O remains in the surface fee sites; at #tot ~ 0.75 ML it 
reaches a saturation value. Such an increase in the work- 
function upon O chemisorption reflects the high elec- 
tronegativity of the adspecies that results in an induced 
inward dipole moment, i.e., with the negative charge at 
the vacuum side of the surface. The saturation at higher 
coverages is then a consequence of the dipole-dipole in- 
teraction giving rise to a depolarization with decreasing 



0-0 distance. 

The initial incorporation of oxygen in the fec/octa and 
hep/octa phases leads to a decrease in the work func- 
tion compared to the saturation value at 1.0 ML with 
O just occupying on-surface fee sites. This decrease is 
considerably smaller for the hep/octa structure, cf. Fig. 
|j; where it should also be noticed that the average bond 
length, O on -Rhi and interlayer spacing, do on -Rhi, be- 
tween chemisorbed O and Rh atoms at the surface are 
practically independent of the occupation of the sites at 
a given coverage, cf. Tables [0] and ILL Moreover, the 
calculated initial-state O Is shifts schematically shown 
in Fig. [|, indicate that the Ols levels of adatoms in the 
hep/octa geometry are ~ 0.5 eV less bound compared 
to corresponding levels of oxygen atoms in the fec/octa 
geometry. This reflects a larger electrostatic repulsion at 
the hep sites, where a somewhat more negatively charged 
O would be adsorbed, which in turn correlates with the 
smaller decrease in the work function for the preferred 
site-switch phase. 

On comparing the values of A$ and AOi s for the 
hep/octa and fec/octa phases with those for the pure 
fee adsorption in Figs. [| and || it can thus be suggested 
that a similarly charged O would sit in hep sites (rather 
than in fee sites) upon octahedral sub-surface occupation 
compared to the oxygens adsorbed in fee sites before O 
penetration__tIence, similarly to the pure on-surface fee 
adsorptionpE3 we suggest that a stronger ionic bonding 
favors the hep sites in the site-switch phases at the onset 
of O penetration (6> tot « 1.0 ML). 

Although we can thus rationalize the higher stability 
of .the hep/octa site-switch phase suggested by Wider et 
al!3 compared to the fec/octa phase, we nevertheless 
stress that the occupation of octahedral sub-surface sites 
is according to our calculations only of a metastable char- 
acter, as the fcc/tetra-I structural combination is con- 
sistently energetically more favorable over the complete 
tested coverage range, cf. Fig. |[ 



C. Continued oxidation: island formation and 
trilayer shift 

We investigate the continued oxidation after sub- 
surface O has been initially incorporated in the coverage 
region 1.25 < 8 to t < 2.0 ML by calculating the aver- 
age adsorption energy at both 1.25 ML (# su b = 0.25 ML) 
and 2.0 ML (9 su h = 1-0 ML) coverages. Again, we find 
that the site combinations hep/octa and fcc/tetra-I give 
the highest adsorption energies, while the occupation of 
tetra-II sites leads to energetically slightly less favorable 
phases, cf. Table I. In the following we will restrict our 
discussion to the two most stable phases, i.e. hep/octa 
and fcc/tetra-I. 
Tables H (6> t , 



Jtot = 1.25ML) and |lll| (0 tot = 2.0ML) 
show that the geometrical changes in response to sub- 
surface penetration of oxygen are significant and it is 
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TABLE II. Calculated structural parameters in A for the 
hcp/octa and fcc/tetra-I phases at (9 tot = 1.25 ML. O atoms 
and all Rh atoms in the two outermost layers were allowed to 
relax. For the interlayer distances, the center of mass of each 
layer is used. Numbers in parenthesis correspond to bulk val- 
ues, which were fixed. O on ,sub-Riii,2 indicate (averaged) bond 
lengths to first and second layer Rh atoms, and AzR ni ,2 and 
Aza on the magnitude of the buckling of the outermost Rh lay- 
ers and of the on-surface O adlayer, respectively. The lateral 
displacements, radially away from the ideal lattice positions, 
are denoted by Arrthi 2 > an d Aro on , respectively. 





hcp/octa 


fcc/tetra-I 


Oon-Rhl 


1.94 


1.95 


S ub-Rhi 


2.15 


2.00 


O sub -Rh 2 


2.21 


1.94 


^Oon-Rhl 


1.14 


1.15 




2.79 


2.78 




2.30 


2.32 




(2.213) 


(2.213) 


Az Rhl 


0.30 


0.23 


Az R h 2 


0.26 


0.30 




0.18 


0.18 


Ar Rhl 


0.06 


0.08 


Ar Rh2 


0.03 


0.00 


Ar 0on 


0.07 


0.07 



TABLE III. Calculated structural parameters in A for the 
hcp/octa, fcc/tetra-I, and fcc/tetra-If au i t phases at 8 to t = 2.0 
ML. O atoms and all Rh atoms in the two outermost layers 
were allowed to relax. Numbers in parenthesis correspond to 
bulk values, which were fixed. Oon.sub-Rhi^ indicate bond 
lengths to first and second layer Rh atoms, and dij interlayer 
spacings. 





hcp/octa 


fcc/tetra-I 


fcc/tetralfauit 


Oon-Rhl 


1.96 


1.97 


1.97 


Oaub-Rhl 


2.10 


2.05 


2.06 


O sub -Rh 2 


2.51 


2.01 


2.02 




1.19 


1.19 


1.19 


di2 


3.36 


3.34 


3.35 


d23 


2.14 


2.17 


2.18 


^34 


(2.213) 


(2.213) 


(2.213) 



therefore not accurate to assume that the Rh(lll) lattice 
will remain essentially undisturbed upon the occupation 
of sub-surface sites. For both hcp/octa and fcc/tetra-I 
phases at 1.25 ML (0.25 ML below the surface), the mean 
outermost substrate interlayer spacing, di2, expands by 
about ~ 22%, and ~ 26%, relative to the corresponding 
value for the bulk-terminated surface, respectively (cf. 
Table [0j) . Prior to the present work, we have calculated 
an increase of di2 from 2.21 A for the clean surface to 
2.37A for the (1 x l)-0/Rh(lll) structure with oxygens 
in fee sites,Q which means that an additional ~ 15 — 
19% expansion is induced by the 0.25 ML of sub-surface 
oxygen. 

For the (2 x 2)-(40 on + O sub )/Rh(lll) phases at 




=1.04 



FIG. 6. Top and side view of the fcc/tetra-I geometry at 
dtot ~ 1.25 ML. Small and large spheres represent oxygen and 
Rh atoms respectively, where those lying in the same plane 
and equivalent under the threefold rotation symmetry have 
the same color. The arrows indicate the direction of reference 
for the atomic in-plane displacements ArR ni 2 , and Aro on of 
Table [n|. Distances are in A. 

1.25 ML, symmetry allows for (in part considerable) 
buckling of the outermost Rh and on-surface O layers, 
i.e. quite distinct local contractions and expansions to- 
gether with lateral shifts (cf. Table |J), which can most 
often be understood as a local expansion oL-the metal 
lattice around the occupied sub-surface site.ES In view 
of keeping the paper within a limited length, we show 
only one example. For instance, in the fcc/tetra-I ge- 
ometry, cf. Fig. |], the first layer Rh atoms above the 
occupied sub-surface site, are pulled out of the surface 
by ~ 0.23 A, while the neighboring Rh atoms in the sec- 
ond substrate layer are pushed inwards by ~ 0.30 A. Also, 
there are local vertical displacements of the O adatoms 
(~ 0.18A), and considerable lateral shifts radially away 
from the ideal lattice positions for both substrate and 
adsorbate atoms. 

The expansion of the first Rh layer distance becomes 
with « 51% even more pronounced at 9 tot = 2.0 ML 
(1ML O below the surface), cf. Table |ffl]. Despite this 
significant distortion of the metal lattice, we find in- 
creased adsorption energies for both the hcp/octa and 
fcc/tetra-I phases compared to the low sub-surface O 
coverage at 0t o t = 

1.25 ML, cf. Fig. g. This reflects 
an attractive interaction between the sub-surface species, 
which implies that the latter have a tendency to form sub- 
surface islands (i.e. nucleate) with a local (1 x 1) period- 
icity. For both the hcp/octa and fcc/tetra-I phase, the 
structure of these islands can be viewed as an O on -Rh- 
O su b trilayer on top of a Rh(lll) substrate, cf. Fig. ||. 
Interestingly, the internal geometry of the strongly bound 
trilayer is almost identical for both phases, cf. Table III 
and Fig. @; it is only the coordination to the underly- 
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fcc/tetral * 




stacking fault 



0.4 0.8 1.2 1.6 
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FIG. 7. Calculated average adsorption energy during the 
registry shift of the O on -Rh-O s ub trilayer along the [211] direc- 
tion over the Rh(lll) substrate. The initial hcp/octa phase 
ends in a fcc/tetra-I configuration with a fault in the stacking 
sequ ence of the surface Rh layer after shifting the trilayer by 
^^372 (an = 3.83 A). 

ing substrate which is different in both cases. In fact, 
we even find that the higher energy hcp/octa configu- 
ration is unstable against a registry shift of this whole 
trilayer along the [211] direction. The calculated average 
adsorption energy along this sliding of the trilayer over 
the Rh(lll) surface is shown in Fig. ^. We have fully 
optimized the structures at the calculated points along 
the barrierless displacement, starting from the hcp/octa 
structure. At the end of the shift by ^-y^3/2 (an = 3.83 
A, Rh lattice constant!]), the on-surface oxygen atoms are 
located in fee and the sub-surface oxygen atoms in tetra- 
I sites, but the surface Rh layer inside the trilayer does 
not continue the fee lattice stacking, but is now located 
in a stacking fault position. The calculated structural 
parameters (interlayer spacings and bond lengths) be- 
tween these two fee / tetra-If au it and fee /tetra-I geometries 
have remained virtually unchanged (cf. Table |l|), and 
the average binding energies differ by only ~ 3 meV/O 
atom, i.e. they are degenerate within our calculational 
uncertainty, which is plausible as the two structures dif- 
fer in fact only by a 60° rotation of the trilayer with 
respect to the underlying substrate. We have similarly 
compared the geometries and average binding energies 
of both, fcc/tetral and fcc/tetralf au it structures at 1.25 
ML and also found no significant difference; we find a 
difference of ~ 1 meV/atom in the calculated binding 
energies. 

The reason behind the ~ 0.2eV/O atom (~ 
0.4eV/unit cell) preference for the fcc/tetra-I structure 
compared to the hcp/octa phase can be found when ana- 
lyzing the coupling of the formed trilayer to the underly- 
ing Rh(lll) substrate. Whereas we consistently find O- 
Rh bondlengths of ~ 2.0 A for both on- and sub-surface 
oxygens in the majority of tested structures, the O su b- 
RI12 bondlength of the sub-surface oxygen to the second 



layer Rh atom in the hcp/octa phase at 6 tot = 2.0 ML is 
significantly larger (2.51 A). This points to a weak cou- 
pling of the trilayer to the Rh(lll) substrate in this ge- 
ometry, which can also be seen in the density plots shown 
in Fig. where the incorporation of sub-surface O in oc- 
tahedral sites is found to induce only very small changes 
on the valence charge of the second layer Rh atoms. 
These changes are considerably larger for the fcc/tetra-I 
configuration, suggesting that the energetic preference of 
this structure is primarily due to an improved coupling 
of the formed trilayer to the underlying substrate. 

Hence, even if sub-surface oxygen was initially incor- 
porated into theanetastable octahedral sites as suggested 
by experiment ,113 it would due to this instability trans- 
form into the fcc/tetra-I configuration upon continued 
oxygen penetration into the Rh(lll) surface. It. is in- 
teresting to notice that on the Ru(0001) surfacetJ'O a 
plausible oxidation pathway has been identified, in which 
after the formation of such trilayers a phase transforma- 
tion into the rutile RuO2(110) structure was obtained for 
local oxygen coverages exceeding 5 ML. The local forma- 
tion of such precursing structures via the agglomeration 
of sub-surface oxygen atoms below the surface that fol- 
low the initial oxygen incorporation, could therefore be a 
more general phenomenon in the oxidation of transition 
metal surfaces. Once the local oxygen coverage exceeds 
a critical value, these precursor configurations will then 
undergo a structural change and actuate the formation 
of bulk oxide phases. In this respect we notice that al- 
ready in the preferred fcc/tetra-I configuration the Rh 
first layer atoms are sixfold coordinated to oxygens and 
the oxygen atoms in the tetrahedral site fourfold coor- 
dinated to metal atoms, i.e. they exhibit identical local 
coordinations as in the most stable, corundum-structured 
RI12O3 bulk oxide. 

IV. SUMMARY 

In conclusion we have presented a density-functional 
theory study addressing the initial penetration of oxy- 
gen into the Rh(lll) surface. Due to the large local 
expansion of the metal lattice induced by the occupa- 
tion of sub-surface sites, O chemisorption on the sur- 
face is initially significantly more favorable. In agree- 
ment with recent experimental findings, we therefore find 
oxygen penetration to occur only after the adsorption of 
almost a full monolayer on the surface. A particularly 
interesting result is that the calculated energy difference 
of ~ 0.3 eV/atom between the pure on-surface and the 
mixed on/sub-surface phases at 6 tot — 1.0 ML is about a 
factor of three smaller than that for Ru(0001). Oxygen 
penetration has at the Ru surface only been reported to 
occur after the (1 x l)-0 phase on the surface has been 
completed. 

The experimentally suggested occupation of octahe- 
dral sites in connection with a site-switch of the on- 
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surface oxygens from fee to hep sites is indeed found to be 
initially possible as a metastable configuration. Increased 
sub-surface O incorporation will then lead to the agglom- 
eration of sub-surface species. The hep/octa phase is 
unstable against a registry shift of the O on -Rh-O su b out- 
ermost layers by which the initial hep/octa phase ends in 
a fcc/tetra-I configuration with a fault in the stacking se- 
quence of the surface Rh layer. In this most stable phase, 
the on- and sub-surface oxygens occupy fee and tetra-I 
sites, respectively, and metal atoms are sixfold coordi- 
nated to oxygens while the sub-surface oxygen atoms in 
tetrahedral sites are fourfold coordinated to metal atoms. 

The similarities between these results and those of a 
simultaneous theoretical study of the oxidation of the 
Ru(0001) surface provide a basis for the interpretation of 
the role of oxygen incorporation, and nucleation as pre- 
cursors of the final bulk oxide structure. Addressing the 
phase transformation to the RI12O3 bulk oxide will thus 
be of considerable future interest, completing the atom- 
istic pathway of oxide formation at the Rh(flf) surface. 



V. ACKNOWLEDGEMENTS 

We thank T. Greber, J. Osterwalder, A. Seitsonen, C. 
Stampfl, F. Illas, and M.E. Grillo for helpful discussions. 



* Present address: Institut fur Chemie, Humboldt Univer- 
sitat zu Berlin, Unter den Linden 6, D-10099 Berlin (Ger- 
many) . 

1 S.H. Oh and J.E. Carpenter, J. Catal. 80, 472 (1983). 

2 G.L. Kellogg, Phys. Rev. Lett. 54, 82 (1985); Surf. Sci. 
171, 359 (1986). 

3 C.H.F. Peden, D.W. Goodman, D.S. Blair, P.J. Berlowitz, 
G.B. Fisher, and S.H. Oh, J. Phys. Chem. 92, 1563 (1988). 

4 A. Bottcher, H. Niehus, S. Schwegmann, H. Over, and G. 
Ertl, J. Phys. Chem. B 101, 11185 (1997). 

5 H. Over, Y.D. Kim, A.P. Seitsonen, S. Wendt, E. Lundgren, 
M. Schmid, P. Varga, A. Morgante, and G. Ertl, Science 
287, 1474 (2000). 

6 Y.D. Kim, H. Over, G. Krabbes, and G. Ertl, Topics in 
Catal. 14, 95 (2001). 

7 G. Comelli, V.R. Dhanak, M. Kiskinova, K.C. Prince, and 
R. Rosei, Surf. Sci. Rep. 32, 165 (1998); and references 
therein. 

8 D. Loffreda, D. Simon, and P. Sautet, J. Chem. Phys. 108, 
6447 (1998). 

9 M.V. Ganduglia-Pirovano and M. Scheffler, Phys. Rev. B 
59, 15533 (1999). 

10 M.V. Ganduglia-Pirovano, M. Scheffler, A. Baraldi, S. 
Lizzit, G. Comelli, G. Paolucci, and R. Rosei, Phys. Rev. 
B 63, 205415 (2001). 

11 P. A. Thiel, J. T. Yates Jr., and W. H. Weinberg, Surf. Sci. 
82, 22 (1979). 



12 M. Rebholz, R. Prins, and N. Kruse, Surf. Sci. 269/270, 
293 (1992). 

13 K.A. Peterlinz and S.J. Sibener, J. Phys. Chem. 99, 2817 
(1995). 

14 K. D. Gibson, J. I. Colonell, and S. J. Sibener, Surf. Sci. 
343, L1155 (1995). 

15 N. M. H. Janssen, A. Schaak, B. E. Nieuwenhuys, and R. 
Imbihl, Surf. Sci. 364, L555 (1996). 

16 J. Wider, T. Greber, E. Wetli, T.J. Kreutz, P. Schwaller, 
and J. Osterwalder, Surf. Sci. 417, 301 (1998); ibid. 432, 
170 (1999). 

17 K.D. Gibson, M. Viste, E.C. Sanchez, and S.J. Sibener, J. 
Chem. Phys. 110, 2757 (1999); ibid. 112, 2470 (2000). 

18 A. Kiejna and B. I. Lundqvist, Phys. Rev. B 63, 085405 
(2001); ibid. Phys. Rev. B 64, 049901 (2001). 

19 K. Reuter, C. Stampfl, M.V. Ganduglia-Pirovano, and M. 
Scheffler, Chem. Phys. Lett. 352, 311 (2002). 

20 K. Reuter, M.V. Ganduglia-Pirovano, C. Stampfl, and M. 
Scheffler, Phys. Rev. B (in press). 

21 J. P. Perdew, J.A. Chevary, S.H. Vosko, K.A. Jackson, M.R. 
Pederson, D.J. Singh, and C. Fiolhais, Phys. Rev. B 46, 
6671 (1992). 

22 P. Blaha, K. Schwarz, and J. Luitz, WIEN97, A Full 
Potential Linearized Augmented Plane Wave Package for 
Calculating Crystal Properties, Karlheinz Schwarz, Techn. 
Universitat Wien, Austria, (1999). ISBN 3-9501031-0-4. 

23 B. Kohler, S. Wilke, M. Scheffler, R. Kouba, and C. 
Ambrosch-Draxl, Comp. Phys. Commun. 94, 31 (1996). 

24 M. Petersen, F. Wagner, L. Hufnagel, M. Scheffler, P. 
Blaha, and K. Schwarz, Comp. Phys. Commun. 126, 294 
(2000). 

25 D.G Castner and GA. Somorjai, Appl. Surf. Sci. 6, 29 
(1980). 

26 M. Todorova, W.X. Li, M.V. Ganduglia-Pirovano, C. 
Stampfl, K. Reuter, and M. Scheffler, Phys. Rev. Lett, (sub- 
mitted) . 

27 The value of « 0.8 eV/O atom was calculated by com- 
paring the average adsorption energy of (2 x 2)-(30 on + 
O su b)/Ru(0001) ordered structures with that of a full ad- 
layer on the surface [(1 x l)-O/Ru(0001), with O on hep 
sites]. The coadsorbate structures considered on Ru(0001) 
are analogous to the fcc/tetral (a)-(b), and the hep/oct (a)- 
(b) of Fig. [lj Calculations were performed within the same 
computational scheme as employed in the present study, 
modelling the metal surface by a six layer slab. Details of 
the used FP-LAPW basis set are those described in Refs. 

Hi 

28 S. Lizzit, A. Baraldi, A. Groso, K. Reuter, M. V. 
Ganduglia-Pirovano, C. Stampfl, M. Scheffler, M. Stich- 
ler, C. Keller, W. Wurth, and D. Menzel, Phys. Rev. B 63, 
205419 (2001). 

29 C. Stampfl, S. Schwegmann, H. Over, M. Scheffler, and G. 
Ertl, Phys. Rev. Lett. 77, 3371 (1996). 

30 A. Bottcher and H. Niehus, J. Chem. Phys. 110, 3186 
(1999). 



8 



(2 x 2) - (40 on + 40 sub )/Rh(lll) tot = 2.0 ML 
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FIG. 8. Left: Contour plots of constant electron density in a [211] plane perpendicular to the (111) surface of 
(2 x 2) — (40 on + 40 su b)/Rh(lll) structures with oxygens in hcp/octa (top) and fcc/tetra-I (bottom) sites. Right: Plot 
of the difference charge density of the same plane, i.e. where the charge density of the clean Rh(lll) surface (with distances 
between Rh atoms as in the chemisorbed system) and those of the isolated oxygens have been subtracted from the electron 
density shown to the left. The distance between contours is 0.3 e/A 3 (left) and 0.06 e/A 3 (right). In both cases, the improved 
bonding of the trilayer to the underlying Rh(lll) substrate can be seen as an increase in the bonding charge density to the 
second layer Rh atoms. 
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